Phase Space Reduction and the Instanton Crossover in (l+l)-Dimensional Turbulence 
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We study (1+1) -dimensional turbulence in the framework of the Martin-Siggia-Rose field theory 
formalism. The analysis is focused on the asymptotic behaviour at the right tail of the probability 
distribution function (pdf) of velocity differences, where shock waves do not contribute. A BRS- 
preserving scheme of phase space reduction, based on the smoothness of the relevant velocity fields, 
leads to an effective theory for a few degrees of freedom. The sum over fluctuations around the 
instanton solution is written as the expectation value of a functional of the time-dependent physical 
fields, which evolve according to a set of Langevin equations. A natural regularization of the 
fluctuation determinant is provided from the fact that the instanton dominates the action for a 
finite time interval. The transition from the turbulent to the instanton dominated regime is related 
to logarithmic corrections to the saddle-point action, manifested on their turn as multiplicative 
power law corrections to the velocity differences pdf. 



I. INTRODUCTION 

As commonly acknowledged through the vast amount 
of numerical investigations of turbulence, intermittency - 
in short, the deviation of pdf 's tails from simple gaussian 
behaviour - turns out to be a phenomenon intimately 
tied to the number of spatial dimensions. The basic phys- 
ical picture behind this observation is that long lived co- 
herent structures, which play roughly the role of local 
reservoirs of conserved quantities, like energy, interact 
more intensely in lower dimensions, perturbing in a rel- 
evant way the cascade process between large and small 
length scales. It is in this sense mostly that Burgers 
model of compressible turbulence in (1+1) dimensions pF 
has been a fruitful ground for theoretical investigation. 
Intermittency is associated there with the existence of 
shock waves, which are described in the inviscid limit, 
v — > 0, by large positive velocity gradients d x u ~ Uq/v 
with support in vanishingly small regions of sizes ~ v/uq. 
Kolmogorov's theory of the inertial range [|2|, which gives 
structure functions S q {r) = (\u(r) — u(0)\ q ) ~ r 2q / 3 , has 
long been known not to match the scaling behaviour 

obtained 



found in the Burgers model, S q (r) ~ ? 
straightforwardly from the fact that the turbulent one- 
dimensional fluid may be depicted as a dilute gas of shock 
waves, interpolated by smooth velocity fields ||,|]]. 

A considerable boost in the understanding of one- 
dimensional turbulence has been achieved in the last few 
years |J4|-jll|, with the help of a number of completely 
different approaches, among them non-perturbative field- 
theory methods. Most efforts have been concentrated on 
the determination of the velocity differences or velocity 
derivatives pdfs. It is now clear that the right and left tail 
of the asymmetric pdfs deserve separate treatments, with 
the left tail results being the subject of some controversy. 

Our aim in this work is to define an effective theory 
for the computation of the right tail pdf, retaining the 
most important degrees of freedom in the analysis. Ini- 



tial attempts on this problem were done by Polyakov jq] 
and Gurarie and Migdal M, who have found leading or- 
der expressions. Our approach closely follows the instan- 
ton method employed by the latter. The computation of 
perturbations around the instanton will be carried out 
here, where, after a natural regularization procedure of 
the fluctuation determinant, subleading logarithm cor- 
rections for the saddle-point action will emerge. The ef- 
fective theory will be defined in such a way as to preserve 
the fundamental BRS symmetry of the Martin-Siggia- 
Rose field theory 0,|l3|, related to general properties 
like renormalization and the existence of a statistical sta- 
tionary state (proved to hold in the class of "purely dis- 
sipative Langevin equations" ) . 

This paper is organized as follows. In sec. II we 
introduce the basic equations to be studied and the 
Martin-Siggia-Rose path-integral formulation. The BRS- 
preserving scheme of phase space reduction is performed, 
and the instanton configuration is defined. Fluctuations 
around the saddle-point fields are then considered in or- 
der to compute the subleading logarithm corrections. In 
sec. Ill we show that, up to first order, viscosity does not 
modify the previous results, in accordance with a well- 
defined inviscid limit. In sec. IV we comment on our 
findings and discuss possible steps of further research. 



II. PHASE SPACE REDUCTION IN THE 
MARTIN-SIGGIA-ROSE FORMALISM 

Let us imagine that the one-dimensional fluid is kept 
in a state of sustained turbulence, by means of the action 
of external random forces. Our starting point is therefore 
the stochastic version of Burgers model, defined by the 
Navier-Stokes equation 

d t u + ud x u = vd x u + / • (1) 

The force / = f(x,t) is taken to be a random gaussian 
variable, with zero mean and two-point correlation func- 



tion 



(f(x,t)f(x',t')} = D(x-x')6(t-t') 



where 



D{x — x') — Dq exp[— (x 



2 /L 2 



(2) 



(3) 



Above, L sets the macroscopic length scale where energy 
pumping is assumed to occur. As L — » oo and the viscos- 
ity vanishes, an intermediate range of wavenumbers forms 
- the inertial range - where scaling behaviour comes into 
play. Arbitrary A-point correlation functions of the ve- 
locity field may be computed in principle through the 
Martin-Siggia-Rose functional HJl3] , 



with 



Z = DiiDuDcDcexp{iS) , 



S= / dxdtu{dtu + ud x u — vd x u) 



(4) 



+ — / dxdx'dtu(x, i)u(x' , t)D(x — x') 
+ / dxdtc(d t c + cd x u + ud x c — vd x c) 



(5) 



Above, c(x, t) and c(x, t) are Grassmann variables. These 
fields enter into the mathematical apparatus as a way to 
account for the non-trivial iacobian defined in the contin- 
uum time formulation JL3 14 1. The Martin-Siggia-Rose 
field theory approach may be used to define a perturba- 
tive expansion involving powers of the convection term, 
where the jacobian's role is just to cancel tadpole dia- 
grams. However, it is important to note that the pertur- 
bative procedure is essentially the Wyld expansion p5| , 
plagued with the well-known insurmountable infrared di- 
vergencies as L — > oo Jl6[ . 

We are interested to study the statistics of the ve- 
locity differences £ = u(p/2,0) — u(—p/2,0) in a non- 
perturbative fashion. The pdf for this observable quan- 
tity may be computed as 



P(0 



1 

2^ 



d\exp(i\t)Z(\) 



(6) 



where 



and 



Z(A)= (exp{-i\[u(p/2,0) - u(-p/2,0)]}> 

DuDuDcDcexp{iS x ) (7) 



(8) 



S X = S- X[u(p/2, 0) - u(-p/2, 0)] . 



We will consider the analytic mapping A — » iX in the next 
computations. Our results, to be obtained in the large 
A limit, will be, therefore, concerned with the right tail 



of the pdf, where shock waves do not contribute. The 
precise mathematical statement is that after the analytic 
mapping the limit A — > oo is related to C ~^ oo in eq. M), 
whenever the right tail pdf decays faster than exp(— c£) 
for any c > 0. To prove it, we write (R) as Z(iX) oc 
J_ d^PiC) cxp(A£). When A — > oo, we will have 



Z(iA)~exp[lnP(Co) + ACo] 



(9) 



Above, we used the usual saddle-point approximation to 
get the resulting asymptotic expression. The parameter 
£o is obtained as the solution of the saddle-point equa- 
tion, 



!'«> 



-P(C)x . 



(10) 



Observe that from the knowledge of the asymptotic be- 
haviour of Z(iX) we are allowed to obtain -P(Co)- All we 
need, thus, is that Co — > oo when A — ► oo, in order to get 
the form of the right tail pdf. 

Taking, for instance, P{£) ~ exp(— c£ Q ) at the right 
tail, it is clear that Z(iX) converges only for a > 1. Fur- 
thermore, we find Co ~ A 1 ^ 0-1 -*, which grows with A, so 
that the large A limit is related in fact to the form of 
the right tail pdf. This is the state of affairs in Burgers 
turbulence, where, as we will see, a = 3. 

The action S is BRS-invariant, according to the fol- 
lowing variations: 



5u = , Su = ec , 8c = , Sc = 



(11) 



Above, e is a constant Grassmann parameter. It is cru- 
cial, in any kind of computational strategy, to work in 
a BRS symmetry preserving scheme, which ensures the 
interpretation of the path integral Z(X) in terms of some 
stochastic dynamical system. This is the fundamental 
guiding principle to be followed, in order to get a re- 
duced form of the Martin-Siggia-Rose functional. Let us 
define the power series 

oo „ 

u{x,t) — y j cr p (t)x p , <J P {t) = I dxx p u(x,t) 
p=o J 

oo „ 

c(x,t) = ^2r] p (t)x p , fj p (t) = / dxx p c(x,t) . (12) 



p=0 



Applying (11) to the above relations, we immediately get 
the "reduccd-BRS" transformations 



S(T. n = 



Sa p = er] p 



Srjp = , 5rjp 



(13) 



The implicit physical motivation in ( |l2| ) is that the rel- 
evant velocity fields involved in the computation of the 
pdf's right tail are smooth. In other words, the above 
transformation to a discrete set of dynamical variables is 
closely related to the Taylor expansion of the stochastic 
force term in the Burgers equation, which is a meaningful 
procedure in the absence of shock waves. 



A phase-space reduction may be implemented by trun- 
cating the series expansions (|12| ) at some arbitrary order. 
Furthermore, if one wants to study the explicit role of vis- 
cous terms, it is necessary to truncate the series beyond 
first order. The simplest choice is to take the expansion 
of u(x,t) up to the second order, which gives 

S\ = dt{a (<j a + a <7i - 2w 2 ) + (ii{(T\ + a 2 + 2ct CT2) 



+0-2(0-2 + 3<7io- 2 ) + i-^-(a 2 - -^000-2 + T5^i + TI^l) 



2 2 3 

2 - u 12^2 + ^ + ^ 

+Vo(i)o + Voo'i - 2^r/ 2 + r)i<r ) + fji(Vi + 2 ??i0i + 277200 
+2i] Q a 2 ) + 172(772 + 37/! 02 + 37720-1) - id(t)Xpai} . 

(14) 

For A = 0, the above action is invariant under the 
reduccd-BRS transformations, as it should be. The 
saddle-point equations, meaningful in the large A limit, 
are 

00 + 0W1 - 2iat 2 + iD a - ijjV2 = , 

(7i + a\ + 2o- o- 2 + i— — -cti = , 

D „ 3D 

0-2 + 3criU2 - 7—2-0-0 + «-7^o- 2 = , 

<$o - °a°i - 2o-io- 2 - 77 77i - 2771772 = , 

0-1 - o-qo-o - 2& x ai - 36-20-2 - 770770 - 2771771 

-3772772 + i\p5(t) = , 

a 2 + 2v& - 26-10-0 - 3o- 2 o-i - 2771770 - 3772771 = , 

7)0 + 7700-1 - 2^2 + 7?io- = , 

771 + 2?7i<ti + 27720-0 + 27700-2 = , 

772 + 37710-2 + 37720-1 = , 

no - fjoo-i - 2771CT2 = , 

m - fj ao - 2771CT1 - 37720-2 = , 

r7 2 + 2^770 - 27/1 o- - 37720-1 = . (15) 

The exact saddle-point equation for u(x,t) shows it 
evolves with "negative viscosity" . Thus, it is necessary, 
in order to solve (|15|), to impose the boundary condition 
ov,(0 + ) = 0. The instanton found by Gurarie and Migdal 
J7J corresponds to take all variables but u\ and u\ equal 
to zero (note that corrections appear at larger orders), 
which yields 

a [°\t) = j t bxR(t), 

af ] \t) = iXpR(t) 2 



(16) 



where 



i?(t) = (l-t/r)- 1 , 



L 2 
D Xp 



(17) 



The fields in the Martin-Siggia-Rose action may be de- 
fined now as perturbations departing from the saddle- 
point solutions, through the replacement 



0-1 

01 



r (o) 



01 

0"! 



(18) 



Substituting these expressions in (jlj) , we obtain a factor- 
ized form for Z(iX), expressed as the product of the dom- 
inant saddle-point contribution and a correction term, 
related to fluctuations around the instanton: 



Z(7A)=exp[75 (A)]Zi(A) , 



(19) 



where 



So (A) 



2D, 



1/2 



3i 



-(Xp) 



3/2 



2 
Z\{\) = I Y\_ D& p Dcr p DfjpDr]pexp[iSi(X)] 
J p=0 



(20) 



and 
S X (X)= 



dt{a (& +0-00-1 



(0) 
o-QO-i 



2va 2 ) 



+0-l(0-l + 0i 



2a 1 (0 Vi 



20-QO-2) + 62(02 + 30102 



, Q (0) \ . .A) /„ 

+ia\ V 2 ) + z — (0, 



2 
L 2 



A 



+770(770 



(0) 
To0i 



21/772 



L 
7?i o-q) 



L*' 



77001 - 

277101 + 2?7i0j u; + 277200 + 277002) 



r (o) 



+772(772 + 377102 + 3?720i + 377201 
+0i O) (0?+ 20 O 2 )}. 



(oh 



(21) 



Observe, as usual, that there are no first order terms 
in the definition of Si (A), since they vanish due to the 
saddle-point equations. The path integral Zi(A) has a 
simple and interesting interpretation. We may write 



Z 1 (X) = (exp[i dta[ '(af 

J —OO 



20002)]) 



(22) 



where the above average is determined by the following 
stochastic (Langevin) equations, 



00 + 0001 + 000-1 - 2t/0 2 = f (t) 

01 + 2oi«tJ 0) + a\ + 20002 = fi(t) 

02 + 30 2 0i + 30102 = h{t) ■ 



(23) 



The time-dependent random forces are correlated as 
(f i (t)f j (t'))=S(t-t')Q ij ,with 



n 



Do 

L 2 



L 2 -1 

2 

-1 3L- 2 



(24) 



The mapping between eqs. ( J20| ) and p^ ) is exact. The 
proof follows from an explicit construction: starting with 
the stochastic system given by eqs. ( p3| ) and (24), we 
find, through the Martin-Siggia-Rose formalism, that the 



expectation value in eq. (|22j) may be precisely written as 
eq. (|20|). This exact mapping is a central point in the 
analysis. 

Taking into account only the linear terms in the 
Langevin equations rt23), in accordance with the theory of 
quadratic fluctuations around the saddle-point, and con- 
sidering the limit of vanishing viscosity, we are able to 
get the exact solution, with initial condition a p {—T) = 
(where it is meant T — * oo) 



Op(*)=exp[-(p + l) f dticftyi)] 
Jo 



x &a/ P (ta)exp[(p + l) / dfeoTte)] 



T 



,(°)/ 



dt 1 f p (t')R(t')^ 



(25) 



A direct application of Wick's theorem shows that inter- 
nal contractions of a\ + la^ai vanish. Only contractions 
among fields defined at different times will produce non- 
vanishing results. We have, thus, expanding the expo- 
nential in (E3), the normal-ordered expression 



v— * * 



i] 



^i(A)=^-/ dt p a[°\tp) 

n 

P =l 



(26) 



There are divergences in the computation of Z\(\) as 
T/t — > oo. Taking the expansion up to second order, it 
is not difficult to get 



Z X (X) ~ ln(T/r) 



(27) 



Reasoning in terms of the cumulant expansion, we find 
that the above logarithmic behaviour clearly implies that 
Z\ scales like a power of T/t. The regularization of this 
divergent result is a subtle point. Note that for times 
t <C — t, we expect fluctuations to be governed by the 
full Martin-Siggia-Rose action describing Burgers turbu- 
lence, since the instanton becomes in this time range es- 
sentially irrelevant. Thus, it is natural to suppose that 
time splits into two regions. We may define in fact a 
time instant to so that for t < t (region I) fluctuations 
of a p are described by pure Burgers turbulence, where as 
the interval to < t < (region II) is dominated by the 
instanton, where it is consistent to apply the theory of 
quadratic fluctuations. The velocity field in region II is 
assumed to be smooth, that is 



cr ~> o\P » a 2 p > .. 



(28) 



As we will see, this condition is satisfied by the solutions 
of the Langevin equations, provided that L 3> p. To find 
to, just take any of the Langevin equations, as the one 
for CTi(i), for instance. Looking at eq. (p3|), we are then 



interested to know at what time the quadratic term a\ 
becomes more relevant than o\o\ . A pragmatical way 
to comp are t he role of these two terms is just to replace 
o\ by \J (cr\) ■ That is, we have to compute 



9(t) 



r (0) 



R(t) 



VM) T^(ri) 



(29) 



The time instant to is defined by the condition g(to) 
Using (E5h we get, in the limit T/t — + oo, 



1. 



2tDo 
3L 2 



R(ty 



(30) 



Assuming now |to/T"| 3> 1, which holds in the asymptotic 
limit of large A, we find 



to 



r2 



(31) 



The basic conclusion drawn from the above arguments is 
that we are enforced to solve ( p3| ) for an arbitrary set of 
initial conditions at t — t Q - Let in this way a p = u p {to) 
represent the initial value of the physical fields in region 
II. In order to match regions I and II, we write 

Z\{\)= / da daida2P(d-o, 0-1,0-2) 

Jc 

x(exp[* / dta i ° ) (al + 2a a 2 )}) . (32) 



The subindex C denotes integration over smooth config- 
urations of the velocity field, and P(oo, o\, 02) is the pdf 
for the set of dynamical variables at time to, which en- 
codes all the information about fluctuations developed 
in region I. According to the picture of Burgers turbu- 
lence as a dilute gas of kinks of the velocity field H , the 
imposition of smoothness in the integration is attained 
when the coordinates x = ±p/2 defined in the evaluation 
of velocity differences are placed in the region between 
shocks. Therefore, a physical interpretation here is to 
regard |in| as the mean time spent at an arbitrary point 
until a shock wave crosses it. 

The solutions of the Langevin eqs. (|23]), which con- 
tain now the explicit dependence on initial conditions, 
are given by 



o P (t) = &(*) + 



R{t ) 

R(t) 



P+i 



(33) 



with 



£ P (i)=i?(i)-^ +1 > / dt'f p (t')R(t')^ . (34) 

It is important to check, as discussed before, the smooth- 
ness condition ( f28[) for the above solutions. We obtain 



M*) a >=«„(*) 2 > + 



R{t ) 



R(t) 



2(p+l) 



<*2> 



n pp i?(i)- 2(p+1) / rfi'i?(i') 2(p+1) 

-'to 

(O , (35) 



i?(i) 

where (ct 2 ) is the average over initial conditions. From 
pure dimensional analysis, we have 

{a 2 p )~(L^/t f . (36) 

Now, since dR(t)/dt > 0, we get, for the integral in (|3£ 

di'i?(0 2(p+1) = / dt'R(t') 2( -P+^R(t'y 



/.. 



(o 



> R(ty 2 / dt'R(t') 2 ( p+ V . (37) 

J to 

Applying this inequality back to ( pq ) and taking into 
account that 



n 



;</' 



n 



p+i,p+i 






L 2 >p 2 



and also 



we get 



R(to) 

[R(t) 



2(p+2) 



< 



fi(«o) 
R(t) 



2(p+l) 



2\„2 



<a p (i) 2 ) » (a p+1 (ty)p- 



(38) 



(39) 



(40) 



which is the smoothness condition fl2Sh stated in a more 
precise way. 

The fluctuating field £(t) is a linear functional of the 
force field, with correlation function 



A lJ (t 1 ,t 2 ) = (Z l {ti)^(t 2 )) = 
= TR(t l )- {i+1) R(t 2 )- (j+1) - 
+F i+j (ti)Q(t2-h)} , 
where 



fi 



(i+j + l) 



'' [F i+j (t 2 )Q(t 1 -t 2 ) 



F p (t)= (p + l)r- J / di'i?(t') (p+2) 
= i?(i) (p+1) - |r/t | (p+1) • 



(42) 



Note that Aij(ti,t 2 ) is a positive definite operator. The 
expression for Zi(X) becomes the product of two factors: 



Zi(X) = Z la (\)Z lb (X) , 



(43) 



where 



/■0 

Zi«(A) = (exp[ ? / d**^? + 2^6)]) 

J to 

= (cx P [- / dt/" dt'& (*)£„(*, *')&(*')]> 

Jto -/to 

^ib(A) = / da daida 2 P(a ,a 1 ,a 2 ) 
Jc 

pO pO 

xexp[/ dt dt'Ji(t)Mij(t,t')J k (t') 

J tn J tr, 



-?-(a 2 + 2a a 2 )}, 



(44) 



with 



B(t,t') = XpR(t) 2 5(t-t') 



J = 2Xp 



R(to) 2 a! 
R(t)R(to)a 



ro 





ii 





f 





i 


n 


n 


) 







and 



M ij (t,t') = [{A- 1 +4B)- 1 ] ij (t,t') 



(45) 



(46) 



We establish now some arguments in order to show 
that Zu has a finite limit as r — ► 0. Two cases have 
to be distinguished here, depending on B being or not a 
positive definite operator. In the second order truncation 
we have been discussing, B is not positive definite. In 
this situation, we were able only to devise a perturbative 
analysis. We write 



M = 



1 



A- 



4B 



[1 - AAB + (AAB) 2 + ...]A . (47) 



Considering the first order approximation to the above 
series, a computer-assisted evaluation (it is necessary to 
examine a large number of terms) shows that both 



pO 

dtx / dt 2 Ji(tl)Aij(tl,t2)Jj(t2) 

to J to 



(48) 



( 41 ) and 



pQ pO pO 

dti / dt 2 / dt$ / dt4 

to J to J to J to 

xJi(ti)A ij (ti,t 2 )B jk (t 2 ,ta)Aki(t 3 ,U)Ji(t4 



(49) 



are finite when r — ► 0. We conjecture that the same 
behaviour holds at any order in perturbation theory. On 
the other hand, if B is a positive definite operator, as 
it happens for the truncation of (O) at first order, it 
is possible to address a proof of the above conjecture. 
At first, we observe that due to the fact that A and B 
are both real, positive definite and symmetric operators, 
inner products satisfy to 



<*, {A- 1 + AB)- 1 ^) < (*, (A + jB- 1 )^) . (50) 

To understand why this is so, just define $ = (A^ 1 + 
AB)- 1 ^. Therefore, 

(tf , (A + -B- 1 )^) = ($, {A- 1 + AB)(A + -B- 1 ) 
x(A- x + 4B)$) = 3($, (A" 1 + 45)$} 
+-{$, J 4" 1 J B" 1 J 4" 1 $) + 16{$,SA J B$) 

= 3($, {A- 1 + 4B)" 1 *) + t^ 1 *, B^A" 1 *} 

+ 16(5$, AB$) > (*, (A- 1 + AB)- 1 ^) . (51) 

This theorem may be applied now to the JMJ term con- 
tained in the definition of Zu,: 

dt / dt'Ji(t)M ij (t,t')J j (t') 

to J to 

< dt dt'J i (t)[A+-B- 1 ] ij (t,t')J j (t') . (52) 
J to Jt 4 

We get in fact a finite limit for the rhs of the above 
expression, in the limit r — ► 0. 

It is interesting to discuss a possible physical interpre- 
tation for the result that Zu,(\) does not depend on A in 
the asymptotic limit. Note that: i) the function Zu,(X) is 
actually defined in terms of the coupling between fluctu- 
ations around the instanton solution and the initial con- 
ditions for the CTi's at time to- If we had fixed Oi(to) = 0, 
for instance, then Zu,(X) would not depend on A, as it 
can be easily derived from eqs. Q) and p5|); ii) recall- 
ing that to is the mean time spent until a shock wave 
crosses an arbitrary point in the one-dimensional space, 
we may regard to hi an equivalent way as the "mem- 
ory time" associated to large and positive fluctuations of 
C (the "right tail" region). In other words, considering 
-f (C2IC1! *2 — ti) - the conditional pdf to observe C2 at 
time t 2 if £1 was observed at time t\ - we have that for 
large values of £2 and for t > \to\ it holds 



P(C 2 |Ci;t)--P(C2) 



(53) 



that is, the conditional pdf becomes the pdf itself at the 
right tail, for t > \to\- Now, as the large C limit is related 
to the large A limit, we find, from i) and ii), that the 
basic physical reason for Zn,(X) not to depend on A (in 
the asymptotic limit) is the absence of memory effects 
from fluctuations which took place in times far enough 
in the past. 

We have, thus, to concentrate our attention on Z\ a . 
Since only quadratic terms are involved in the definition 
of Z\ a , we find 



Zi = exp[--Trln(l + 4^B)] 



(54) 



the first non-zero contribution gives, in the limit r — > 0, 
and after a lenghty computation, 

Z la = exp[4Tr(ABAB)] - exp[cln(|t /r|)] , (55) 

with c = 224/45 ~ 5. Since r ~ A~ 1/2 , we get 

P(C)~C c exp(-C 3 )- (56) 

This is our central result. A careful analysis of the in- 
stanton crossover has led ultimately to the power law 
correction shown in eq. (56), for the asymptotic form of 
the velocity differences pdf at the right tail. 



III. VISCOSITY EFFECTS 

In order to study viscosity effects, all we have to do, 
as it follows from eq. (|23|) , is to perform the replacement 
/o -* /o + 2^cr 2 , that is, 



fo{t) -» /o(*) + 2 v R{t)- A / dt'f 2 (t')R(t'Y + 2ua 2 

J t 



(57) 



We have now, 



(/o(«)/o(f )> = 

= Q 00 S(t - t') + 2M 02 R(t') 3 R(t)- 3 + 0{v 2 ) 

(/o(W)> = 

= Q 02 5(t - t') + 2M 22 R(t') 3 R(t)- 3 0(t - t') . (58) 

The above definitions may be used to compute correc- 
tions to the correlation functions Aij(t,t'). We get, up 
to first order, 



Tr ln(l + AAB) = 4Tr(AB) 



via 

L 2 



(59) 



Expanding the logarithm as ln(l + AAB) — AAB — 
8(ABAB) + ..., and taking into account that Tr(AB) = 0, 



We find, thus, that viscosity does not lead to logarithm 
corrections, and that the limit v — > is a well-defined 
one, at least to first non-trivial order in perturbation the- 
ory. 



IV. CONCLUSIONS 

We studied the instanton crossover problem in the 
field-theoretical approach to Burgers turbulence. With 
to = — (L 2 /Do) 1 / 3 , this crossover is characterized by 
the transition from fully developed turbulence, at times 
t < to, to the instanton dominated regime, in the in- 
terval to < t < 0, where quadratic fluctuations may be 
integrated out. The essential outcome of our analysis 
is that pdf's tails get multiplicative power law correc- 
tions to the leading, saddle-point expressions. Such phe- 
nomenon, which is hard to detect numerically, is likely 
to occur in other intermittent systems. 



We note that some questions are left open, which re- 
quire even more labourious computations. As already 
observed, the operator B, introduced in eq. Ma), may be 
positive definite or not, according to the order where the 
truncation of ( |l2| ) is done. Since the sign of the parameter 
c (defined in eq. (|55|)) is positive if B is positive definite, 
but the converse is not necessarily true, it is clear that as 
higher order truncations are considered, we expect that 
i) B becomes positive definite or not, with c — ► c > 0, or 
ii) B becomes non-positive, with c — > c < 0. We do not 
have enough information yet to decide which one is the 
correct choice; further work is necessary on this matter. 
Also, it is worth mentioning that power law corrections 
may be model-dependent, as strongly suggested by the 
leading order saddle-point results themselves, related to 
the specific form of the force-force correlation function. 

A natural question is if analogous results hold for 
the left tail pdf, where the Martin-Siggia-Rose formal- 
ism may as well be employed |9|. The techniques dis- 
cussed here are adequate to study this problem, a sub- 
ject deserved for future investigations. Even though the 
path integral approach is direct and appealing, it relies 
on heavy perturbative computations performed on trun- 
cated series. It would be interesting, and perhaps use- 
ful, to know how the power law corrections could be de- 
rived from some alternative analysis of randomly forced 
Burgers turbulence (after completion of this work we be- 
came aware of rcf . |17| , where power law corrections are 
rigourously obtained in the case of decaying Burgers tur- 
bulence) . 
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